Generalized Gompertz-power series distributions

arXiv:1508.07634vl [stat.ME] 30 Aug 2015

Generalized Gompertz-power series distributions

S. Tahmasebi^ , A. A. Jafari^’*

'^Department of Statistics, Persian Gulf University, Bushehr, Iran
^Department of Statistics, Yazd University, Yazd, Iran

Abstract

In this paper, we introduce the generalized Gompertz-power series class of distribu¬
tions which is obtained by compounding generalized Gompertz and power series distribu-
tions. This comp oundi ng procedure follows same way that was previously carried out by

Silva et al

(2013) and

Barreto-Souza et al.

(|201lll in introducing the compound class of

extended Weibull-power series distribution and the Weibull-geometric distribution, respec¬
tively. This distribution contains several lifetime models such as generalized Gompertz, gen¬
eralized Gompertz-geometric, generalized Gompertz-poisson, generalized Gompertz-binomial
distribution, and generalized Gompertz-logarithmic distribution as special cases. The haz¬
ard rate function of the new class of distributions can be increasing, decreasing and bathtub¬
shaped. We obtain several properties of this distribution such as its probability density
function, Shannon entropy, its mean residual life and failure rate functions, quantiles and
moments. The maximum likelihood estimation procedure via a EM-algorithm is presented,
and sub-models of the distribution are studied in details.

Keywords: EM algorithm, Generalized Gompertz distribution, Maximum likelihood es¬
timation, Power series distributions.

1 Introduction

The exponential dist ribution i s com monly used in many applied problems, particularly in life¬

time data analysis (jLawlessl .

20031 ). A generalization of this distribution is the Gompertz

distribution. It is a lifetime distribution and is often applied to describe the distribution of
adult life spans by actuaries and demographers. The Gompertz distribution is considered for
the analysis of sur vival in some sciences such as biology, gerontology, computer, and marketing
science. Recently, iGuota and Kundul (|l999l i defined the generalized exponential distribution

* Corresponding: aajafari@yazd.ac.ir

and in similar manner. lEl-Goharv et alJ (|2013l i introduced the generalized Gompertz (GG) dis¬
tribution. A random variable X is said to have a GG distribution denoted by GG(a,/ 3 , 7 ), if
its cumulative distribution function (cdf) is

G(x} = [1 — a,/3 > 0, 7 > 0; x > 0.

and the probability density function (pdf) is

( 1 . 1 )

( 1 . 2 )

The GG distribution is a flexible distribution that can be skewed to the right and to the
left, and the well-known distributions are special cases of this distribution; the generalized
exponential proposed by

Gupta and Kundul (jl9^ 1 when 7 —)• O'*', the Gompertz distribution

when 0 = 1, and the exponential distribution when a = 1 and 7 —?■ O'*'.

In this paper, we compound the generalized Gompertz and power series distributions,
and introduce a new class of distribution. This procedure follows similar way that was previ-
ously carried out by some aut hors: The exponential-power series distribution is introduced by

Ghah

^nnJ

candi and Ganialil (1200911 w

Adamidis and Loukas

(Tahmasbi and Rezaei.

by

2008

rich is concluded the e xponential

1998|), exponential-Poisson (|Kua . 120071 1. and exponential-logarithmic

-geometric (jAdamidis et al.

distri

rutions; the Weibull-power series distributions is introduced

Morais and Barreto-Souzal (j201lh and is a generalization of the exponential- power series dis¬

tribut ion; the generalized exponential-power series di stribution is in t roduced by

2012l i which is concluded the Poisson-expone ntial (iGancho et al

201 ih complementary ex ponential-geometric (jLouzada et al

exponentia’

-power series (IFlores et al

tributions (jMahmoudi and .Tafari

201 4l i.

2011

Mahmoudi and Jafari

Louzada-Neto et al.

201 ih . and the complementary

2013 ) distributions; linear failure rate-power series dis-

The remainder of our paper is organized as follows: In Section [2l we give the probability
density and failure rate functions of the new distribution. Some properties such as quantiles,
moments, order statistics. Shannon entropy and mean residual life are given in Section O In
Section [H we consider four special cases of this new distribution. We discuss estimation by
maximum likelihood and provide an expression for Fisher’s information matrix in Section [5j A
simulation study is performed in Section [H An application is given in the Section [71

2 The generalized Gompertz-power series model

A discrete random variable, is a member of power series distributions (truncated at zero) if
its probability mass function is given by

Pn = P{N = n) =

anO^

c{e)

, n = 1,2 ,...,

( 2 . 1 )

where a„ > 0 depends only on n, C{9) = and 6 E (0, s) (s can be oo) is such that

C{9) is finite. Table [T] summarizes some particular cases of the truncated (at zero) power series

distributions (geometric, Poisson,
series distribution can be found in

ogarit

i mic a nd binomial). Detailed properties of power

Noack ( 1950 b Here, C'{9), C''{9) and C'''{9) denote the

first, second and third derivatives of C{9) with respect to 9, respectively.

Table 1: Useful quantities for some power series distributions.

Distribution

an 0(9) C'(9) C"(9) 0'''(9) s

Geometric

Poisson

Logarithmic

Binomial

1 9(1-9)-^ (1 - 6»)-2 2(1 - 6»)-3 6(1-6 »)-^ 1

n\~^ — 1 e® oo

n-^ -log(l-0) (1-9)-^ (1-0)”^ 2(1-9)-^ 1

('in\ 1 a\m 1 m m(m-l) m(m-l)(fc-2)

\nJ G “T^

We define generalized Gompertz-Power Series (GGPS) class of distributions denoted as
GGPS(a, /3, 7 , 9) with cdf

T(x) f^ ^r.{eG{x)r C{9G{x)) C{9n ^ ^ ^

n=l

C{9)

C{9) C{9)

( 2 . 2 )

where t = I — e t"

. The pdf of this distribution is given by

(2.3)

This class of distribution is obtained by compounding the Gompertz distribution and power
series class of distributions as follows. Let be a random variable denoting the number of
failure causes which it is a member of power series distributions (truncated at zero). For
given N, let Xi,X 2 , ■ ■ ■ ,Xn' be a independent random sample of size N from a GG{a, 13,'y)
distribution. Let = maxi<j<jv A"*. Then, the conditional cdf of | A = n is given by

_ [1 _ p-f

which has GG(na,/3,7) distribution. Hence, we obtain

P(A(7 v) < x,N = n) =

ani9G{x))^ ^ an9^ _ £(e7^_i) „

0(9) Ci9)'- ‘ ■

Therefore, the marginal cdf of has GGPS distribution. This class of distributions can

be applied to reliability problems. Therefore, some of its properties are investigated in the
following.

Proposition 1. The pdf of GGPS class can be expressed as infinite linear combination of pdf
of order distribution, i.e. it can he written as

OO

f{x) = '^Png{n){x;na,P,-f), (2.4)

71=1

where g(„)(x; no,/3, 7 ) is the pdf of GG{na,

— §. — V)

Proof. Gonsider t = 1 — e t" ^ h So

^ /7

= = '^Png{n){x;na,ld,j).

71 = 1 ^ 71=1

□

Proposition 2. The limiting distribution 0 / GGPS(a,/3, 7 , 0) when 0 —)• O'*" is

lim F{x) = [1 — e i)jca

6 ^- 0 +

which is a GG distribution with parameters ca, fd, and 7 , where c = min{n G N : o„ > 0}.

Proof. Gonsider t = 1 — e

. So

lim F(x)
9 ^ 0 +

C(At"'

el^+~C{9)

lim

A->-0+

OO

71=1

OO

E anO^

71=1

= lim

0 ^. 0 +

OO

aci“+ E

n=c+l

OO

flc + E CLnO"-~'^

n=c+l

= F

□

Proposition 3. The limiting distribution 0 / GGPS(a, /?, 7 , 0) when 7 —?■ O"*" is

lim F(x)

7 ^ 0 +

C( 0(1 -

W)

i.e. the cdf of the g e neral ized exponential-power series class of distribution introduced by
Mahmoudi and .Jafari f20lA ).

Proof. When 7 ^ 0 ^, the generalized Gompertz distribution becomes to generalized exponen¬
tial distribution. Therefore, proof is obvious. □

Proposition 4. The hazard rate function of the GGPS class of distributions is

h{x)

eal3e'^^{l-t)t^-^C'{et°^)

C{e) -

(2.5)

ILi.

where t = 1 — e r

(eT'^-l)

Proof. Using ( 12 . 21 ) . (|2.3I) and definition of hazard rate function as h{x) = f{x)/{l — F{x)), the
proof is obvious. □

Proposition 5. For the pdf in (12.3p . we have

lim /(x)
X—> 0 +

C (0)9(3
C(9)

0 < a < 1
0 = 1
a > 1,

lim f{x) = 0 .

X^OO

Proof. The proof is a forward calculation using the following limits

lim

a:->- 0 +

'00 0 < a < 1

< 1 0 = 1
_ 0 o > 1 ,

lim F = 0 , lim t = 1 .
a;->- 0 + x^oo

□

Proposition 6. For the hazard rate function in (|2.5p . we have

lim h{x)

3;^0+

00

c'(o)ep

0 ( 9 )

0 < o < 1
0 = 1
o > 1 ,

lim h{x)

x—>-oo

00 7 > 0
/3 7 —)• 0

Proof. Since lim 3 ,^o+(l “ F{x)) = 1, we have lim 2 ,_^Q+ h{x) = lim 2 ,_^Q+ f{x).
For lima;_>.oo h{x), the proof is satished using the limits

lim c'{er) = c'(e), lim = i,

X—)-oo x^oo

lim = lim ^ f 7 > 0

x^co C(9) — C(Otr x^oo 0l3aC'{6)e^^{l — t) y 9aC'(9) 7 0.

□

As a example, we consider C {9) = 6 + 9“^^. The plots of pdf and hazard rate function of
GGPS for parameters /3 = 1,7 = .01,0 = 1.0, and o = 0.1,0.5,1.0, 2.0 are given in Figured]

This pdf is bimodal when o = 2.0, and the values of modes are 0.7 and 3.51.

Density

Hazard

Figure 1: Plots of pdf and hazard rate functions of GGPS with C {9) = 6 -\- 6“^^.

3 Statistical properties

In this section, some properties of GGPS distribution such as quantiles, moments, order statis¬
tics, Shannon entropy and mean residual life are obtained.

3.1 Quantiles and Moments

The quantile q of GGPS is given by

X

^ ^ 9 ^

0 < O' < 1,

where G~^{y) = ^log[l — ) ] and C'“^(.) is the inverse function of C{.). This result

helps in simulating data from the GGPS distribution with generating uniform distribution data.

For checking the consistency of the simulating data set form GGPS distribution, the his¬
togram for a generated data set with size 100 and the exact pdf of GGPS with C (9) = 9 + 9'^^,
and parameters a = 2, /3 = 1, 7 = 0.01, 9 = 1.0, are displayed in Figure [2] (left). Also, the
empirical cdf and the exact cdf are given in Figure [2] (right).

Consider X ~ GGPS(a,/3, 7 ,0). Then the Laplace transform of the GGPS class can be
expressed as

L{s) = E{e-^^) = Y^P{N = n)Ln{s),

n=l

where Ln{s) is the Laplace transform of GG(na,/ 3 , 7 ) distribution given as

r+oo

LJs) =

f

= na/J

I

E

j=0

(3.1)

Histogram

Empirical Distribution

Figure 2: The histogram of a generated data set with the pdf (left) and the empirical cdf with
cdf (right) of GGPS distribution.

na — 1

j

{-lye

na/3j; ("Y‘) /

{'1-s)x^=y(j+l)e^-

dx

T{k + l)

nafd^Yl

na

-l\ (-l)^'+^ef(^'+')[^(j + l)]

J

T{k + l)(s — 7 — 'jk)

s > 7-

(3.2)

j=0 k=0

Now, we obtain the moment generating function of GGPS.

OO

Mx{t) = E{y^) = Y,P{N = n)Ln{-t)

n=l

n=l ' ' k=0j=0

- n(’^7')(-l)^+"+7^^^+'7f(j + 1))*

T{k + l)(t + 7 + 'jk)

S.I

- “ iV(^“-i)(-l)Wie7(^'+i)(^(j + 1))^

fc=0 j=0

(3.3)

r(A: + l)(t + 7 + 7A;)

where is a random variable from the power series family with the probability mass function
in (EU) and E]y[U] is expectation of U with respect to random variable N.

We can use Mx{t) to obtain the non-central moments, = ElX'"]. But from the direct
calculation, we have

n=l ' ' k=0 j=0

al3Ex[^'^

k=0 j=0

r(k -|- 1)(7 + jky~^^

^(7Va-l)(_^)j+fc+r+lgf (i+l)(^(^- ^ + 1)

r(k + i)(j + jky+^ ^

(3.4)

Proposition 7. For non-central moment function in\3.4[ we have

lim fij. = E[Y^],
e^o+

where Y has GG(ca, (3, 7 ) and c = min{n G N : > 0}.

Proof. If y has GG(ca,/ 3 , 7 ), then

^ ^ hh r{k + i){^ + jky+^

Therefore,

E ane^E[Yy

lim fir = lim ^ -

0 ^ 0 + 0 ^ 0 + ^

n=l

= lim
e^-o+

00

a,E[Yy+ E

n=cH-l

00

«c + I] anO^-^
n=c+l

= E[Yy.

□

3.2 Order statistic

Let Xi,X 2 , ..., Xn be an independent random sample of size n from GGPS(a, (3, 7 , 9). Then,
the pdf of the ith order statistic, say Xi-n, is given by

f. (x) = _ — _ V ~Ml \ n-i

h:n[ ) ^i_i)i^n-i)r^c{e) ^ ^ c{9) ^ ’

_ _§.^Qyx_-^\

where / is the pdf given in ( 12 .31) and t = 1 — e r 7 Also, the cdf of Xi^n is given by

F ui G (-i)n 7 ‘) , c(t“)

' (i-l)!(n-i)!^ k + i+l 7(9)

k-\-i

An analytical expression for rth non-central moment of order statistics Xi-n is obtained as

/t _ 1 \ /r)\ r+°°

fc-n-l-j-l / ^ ^\ '''

k^—Tt — 1~\~ 1

= r

E

(_l)fc-n+*-l _ 1

k=n—i+l

[C{9yy \n — ij \k

:)r

x^-yc{d) - c{ee)fdx

where S{x) = 1 — E{x) is the survival function of GGPS distribution.

3.3 Shannon entropy and mean residual life

If X is a no ne-negativ e cont inuous random variable with pdf /, then Shannon’s entropy of X

is defined by

Sha,Tinonl ( 19481 ) as

^- 1-00

H{f) = E[- log f{X)] = - f{x) log{f{x))dx,

Jo

and this is usually referred to as the continuous entropy (or differential entropy). An explicit
expression of Shannon entropy for GGPS distribution is obtained as

Hif) = E{- log[^e^^(e-f- e'f(^(1 - g-f

— (a — l)ili[log(l — — ili[log(C''^0(1 —

= - “ 71^1 + -Mxil) - - - {a - 1)Y^ P{N = n) f log{t)dt

E{d) 7 7 ^ Jo

oo „!

— ^2 P{X = n) nvP'~^ \og{C'{0u))du

n=i Jo

= -l<>gl^|-Wi + -Afe(7)-- + 7^EN|7]-£„[A(Ar,())], (3,5)

0 ( 0 ) 7 7 a A

where A{N, 0) = Nu^~^ \og{C'{6u))du, A is a random variable from the power series family

with the probability mass function in (12.ip . and Ex\U] is expectation of U with respect to
random variable N. In reliability theory and survival analysis, X usually denotes a duration
such as the lifetime. The residual lifetime of the system when it is still operating at time s, is
Xs = X — s\X>s which has pdf

f{x]s)

fix) eg{x)C'{eGix))
l-F(s) C( 0 )-C( 0 G(s))’ -

Also, the mean residual lifetime of Xg is given by

mis) = E[X - s\X > s]

C°°ix - s)fix)dx
1-Fis)
xfix)dx
l-Eis) ""
Cie)EN[Zis,N)]
C(0)-C(0[l-e“^^"^"“^^]“)

where Z(s,n) = xg^n)(x','na,/3,j)dx, and g(„)(x; na,/1, 7 ) is the pdf of GG(na,/3, 7 ).

4 Special cases of GGPS distribution

In this Section, we consider four special cases of the GGPS distribution. To simplify, we
consider t = 1 — x > 0, and Aj =

4.1 Generalized Gompertz-geometric distribution

The geometric distribution (truncated at zero) is a special case of power series distributions
with On = 1 and C{0) = (0 < 9 < 1). The pdf and hazard rate function of generalized

Gompertz-geometric (GGG) distribution is given respectively by

(l-0)a/3eT'"=(l-f)r-i

Remark 4.1. Consider

fix) =

h{x) =

fuix) =

{9t^ - 1)2

(l-0)a/3eT'^(l-t)t“-i
(1 - 6 »t")(l - t“) '

X > 0,

X > 0.

X > 0,

(4.1)

(4.2)

(4.3)

((1 - 6 »*)t“ - 1)2 ’

where 9* = 1 — 9. Then fuix) is pdf for all 9* > 0 (see
that when a = 1 and 7 —0+, the pdf of extended exponential geometric (EEC) distribution

Marshall and Olkin.

1991). Note

l Adamidis et al\. 12003) is concluded from (14.311 . The EEC hazard function is monotonically

increasing for 9* > 1; decreasing for 0 < 0* < 1 and constant for 9* = 1.

Remark 4.2. If a = 9* = 1, then the pdf in (14.3p becomes the pdf of Gompertz distribution.
Note that the hazard rate function of Gompertz distribution is h{x) = jde"^^ which is increasing.

The plots of pdf and hazard rate function of GGG for different values of a, /3, 7 and 9*
are given in Figure [3l

Theorem 4.1. Consider the GGG hazard function in (SSI). Then, for a > 1, the hazard
function is increasing and for 0 < a < 1 , is decreasing and bathtub shaped.

Proof. See Appendix A.l.

The first and second non-central moments of GGG are given by

□

Em - - 0 E n."- E E +1))‘

k=0j=0

n=l

T{k + 1)(7 -I- 7A:)2

E{X^) = 2a/3{l -9)Y^ n9^-^ EE

k=0j=0

n=l

7L.(_l)i+fe+3efO+i)(|(j + l))fc

r(A:-I- 1)(7 + 7 ^)^

Density

Density

Density

Figure 3: Plots of pdf and hazard rate function of GGG for different values a, /?, 7 and 9*.

4.2 Generalized Gompertz-Poisson distribution

The Poisson distribution (truncated at zero) is a special case of power series distributions
with On = ^ and C{9) = — 1 {9 > Q). The pdf and hazard rate function of generalized

Gompertz-Poisson (GGP) distribution are given respectively by

f{x)

h{x)

9al3e^^-\l-t)e-^e^^\ x>0

0a/3eT'*(l

g6» _ g6>t“ ’ X > U.

(4.4)

(4.5)

Theorem 4.2. Consider the GGP hazard function in ra- Then, for a > 1, the hazard
function is increasing and for 0 < a < 1, is decreasing and bathtub shaped.

Proof. See Appendix A. 2 .

□

The hrst and second non-central moments of GGP can be computed as

- - 7l,-(-l)^'+^ef(^'+^)(f(j + 1 ))'=

^ (n- 1 )! ^ ^

=1 ^ ^ k= 0 j =0

T{k -I- 1)(7 -I-

r 7, .7 (i +1))*

- 1 (™- 1 )! P{k + l){-i+ -ikf

Density

Density

Density

Figure 4: Plots of pdf and hazard rate function of GGP for different values a, /3, 7 and 6.

The plots of pdf and hazard rate function of GGP for different values of a, /3, 7 and 9 are given
in Figured!

4.3 Generalized Gompertz-binomial distribution

The binomial distribution (truncated at zero) is a special case of power series distributions with
an = (™) and C{6) = (0 + 1)™ — 1 (0 > 0), where m {n < m) is the number of replicas. The
pdf and hazard rate function of generalized Gompertz-binomial (GGB) distribution are given
respectively by

fix)

(g+oi _j_ 1 \m—1

= x > 0 ,

(4.6)

h{x)

m0a/3eT'^(l - + 1 )”^-^

“ {9 + 1 )™ - ( 0 t" -Fl)™ , x> .

(4.7)

The plots of pdf and hazard rate function of GGB for m = 4, and different values of a,
(3, 7 and 0 are given in Figure O We can find that the GGP distribution can be obtained as
limiting of GGB distribution if mO —?■ A > 0, when m —?• 00 .

Theorem 4.3.

Consider the GGB hazard function in

Then, for a > 1, the hazard

function is increasing and for 0 < a < 1, is decreasing and bathtub shaped.

Density

Density

Density

Figure 5: Plots of pdf and hazard rate function of GGB for m = 5, and different values a, (3,
7 and 9.

Proof. The proof is omitted, since 9 > 0 and therefore the proof is similar to the proof of
Theorem 14.21 □

The first and second non-central moments of GGB are given by

E(X)

E(X^)

fm] + 1 ))"

(0 + 1 )”*-!^^ r(k + i)(j + jkf

2a/3 ^ fm\ ^ ^ + 1))"

^\nj r(k + l)(j + jk)3

4.4 Generalized Gompertz-logarithmic distribution

The logarithmic distribution (truncated at zero) is also a special case of power series distribu¬
tions with On = ^ and C(9) = — log(l — 9) (0 < 9 < 1). The pdf and hazard rate function of
generalized Gompertz-logarithmic (GGL) distribution are given respectively by

f(x)

h{x)

9ape^^{l-t)t^-^

— 1) log(l — 9) ’

9aPe^^{l-t)t^-^
{9t^ - l)log(^^)’

X > 0,

X > 0.

(4.8)

(4.9)

Density

Density

Density

Figure 6 : Plots of pdf and hazard rate function of GGL for different values a, /?, 7 and 9.

The plots of pdf and hazard rate function of GGL for different values of a, j3, 7 and 9 are
given in Figured

Theorem 4.4. Consider the GGL hazard function in (JfW. Then, for a > 1, the hazard
function is increasing and for 0 < a < 1 , is decreasing and bathtub shaped.

Proof. The proof is omitted, since 0 < 0 < 1 and therefore the proof is similar to the proof of
Theorem 1. □

The hrst and second non-central moments of GGL are

EiX) =

afd

— logfl — 9) ^

J n=l fc=0 3=0

r{k + 1)(7 -I- 7^)2

-log(l-0)^ F(/i; +1)(7 + 7A:)3

5 Estimation and inference

In this section, we will derive the maximum likelihood estimators (MLE) of the unknown
parameters 0 = {a, (3,^,9)'^ of the GGPS(a,/3, 7 , 0). Also, asymptotic conhdence intervals of

these parameters will be derived based on the Fisher information. At the end, we proposed an
Expectation-Maximization (EM) algorithm for estimating the parameters.

5.1 MLE for parameters

Let Xi,...,Xn be an independent random sample, with observed values xi,...,Xn from
GGPS(q;,/3, 7 , 0) and 0 = (a,/ 3 , 7 , 0 )^ be a parameter vector. The log-likelihood function
is given by

n n

ln = ln{®]x) = nlog( 6 ')-t nlog(Q;/3)-I-n 7 x-F ^log(l - L)-h (a - 1 ) ^log(ti)

i=l i=l

n

+ ^log(C'( 0 tf))-nlog(C( 0 )),
i=l

where = 1 — e Therefore, the score function is given by [7(0; x) = (^,

where

din
da

din

dp

din

d'j

din

d9

The MLE of 0, say 0, is obtained by solving the nonlinear system U{Q;x) = 0. We cannot
get an explicit form for this nonlinear system of equations and they can be calculated by using
a numerical method, like the Newton method or the bisection method.

Eor each element of the power series distributions (geometric, Poisson, logarithmic and
binomial), we have the following theorems for the MLE of parameters:

Theorem 5.1. Let gi(a;/?, 7 , 0, x) denote the funetion on RHS of the expression in dSU),
where p, 7 and 6 are the true values of the parameters. Then, for a given /3 > 0, 7 > 0 and
9 > 0, the roots 0 /gi(a,/3; 7 , 0, a;) = 0, lies in the interval

n n

c'(e) + ^ i=i *=i

n

a

+ '^log{ti) +Y^

i=l

i=l

9tflog{U)C"{9tf)
C'{9tf)

1 ^

_

n

n dti n nditf)

tj.

2 = 1

2=1

2=1 2=1

Xie

'yxi^

n dti n

+E ^

n

+ E

i=i ‘ i=i

tfC"{9tf) nC'{9)

C'{9tf

i=l

cm

C{9)

cm

(5.1)

(5.2)

(5.3)

(5.4)

Proof. See Appendix B.l.

□

Theorem 5.2. Let g 2 (/ 3 ; a, 7 , 0, x) denote the function on RHS of the expression in ()5.3p .
where a, 7 and 9 are the true values of the parameters. Then, the equation g 2 (/ 3 ; a, 7 , 9,x) = 0
has at least one root.

Proof. See Appendix B.2. □

Theorem 5.3. Aei g 3 ( 0 ; a,/3, 7 , x) denote the function on RHS of the expression in (|5.4p and
X = n~^ where a, (5 and 7 are the true values of the parameters.

a) The equation g 3 ( 0 ; a,/3, 7 , x) = 0 has at least one root for all GGG, GGP and GGL distri¬
butions if Yl'i=i > i •

b) If g 3 (p; a, ft, 7 , x) = where p = and p G (0,1) then the equation gz{p', a, /3, 7 , x) = 0

has at least one root for GGB distribution if X]r=i •

Proof. See Appendix B.3. □

Now, we derive asymptotic confidence intervals for the parameters of GGPS distribution

It is well-known that under regularity conditions (see

Casella and Bergen.

2001

, Section 10),

the asymptotic distribution of ^/n(& — 0) is multivariate normal with mean 0 and variance-
covariance matrix J“^(0), where Jn(0) = lim„_).oo -ln(0)) and 1^(0) is the 4x4 observed
information matrix, i.e.

I O'a

I a 13

laO

In ( 0 ) = -

IjBo

Ipfi

1-10

ip-1

R-f

Ipe

I-yO

_ loa

lep

IQ"!

lee _

whose elements are given in Appendix C. Therefore, an 100(1—r/) asymptotic confidence interval
for each parameter, 0 ^, is given by

AC R — (0r — ^rj/2 ^ ; 0r 4 “ I,

where Irr is the (r, r) diagonal element of /„^(0) for r = 1,2,3,4 and Z ^/2 is the quantile ^
of the standard normal distribution.

5.2 EM-algorithm

The traditional methods to obtain the MLE of parameters are numerical methods by solving the
equations (|5.in - (l5.4l) . and sensitive to the initial values. Therefore, we develop an Expectation-
Maximization (EM) algorithm to obtain the MLE of parameters. It is an iterative metho d,

and is a very powerful tool in handling the incomplete data problem (jPempster et al

19771).

We define a hypothetical complete-data distribution with a joint pdf in the form

5 ((x,2;;0) = S

where t = 1 — and a,/3, 7 , 0 > 0, x > 0 and z G N. Suppose

is the current estimate (in the rth iteration) of ©. Then, the E-step of an EM cycle
requires the expectation of {Z\X;@^'^'^). The pdf of Z given X = x is given by

g{z\x)

and since

c'{e) + ec"{e) = ^ a^ze^-^ + ^

2=1 2 = 1

'^z^azO'' S

Z = 1

the expected value of Z\X = x is obtained as

E{Z\X = x)

(5.5)

By using the MLE over 0, with the missing Z’s replaced by their conditional expecta¬
tions given above, the M-step of EM cycle is completed. Therefore, the log-likelihood for the
complete-data is

n n

l*n{y,®) OC ^ Zi log(0) -I- n log(a/3) -|- nyx -|- ^ log(l — ti)

i=\ i=l

n

+ '^{zia - l)log(fj) - nlog(C( 6 ')), (5.6)

i=l

where y = (x;^;), x = (xi,...,x,i) and 2 ; = (zi,..., Zn)- On differentiation of (|5.6n with
respect to parameters a, 13, 7 and 6, we obtain the components of the score function, U{y, 0) =

V «

da ’ 8/3 ’’ d'Y ’’ 89 ■' ’

as

di:

_ g 7

2 = 1

dlZ n 1

X =

i) '^{zia - 1)

dl*n P

dl*

'-"‘n

86

E

2=1

2=1

c'{e)

6 ”C( 0 )'

2=1

7

i(e7-. - 1)

n n _ Ij -I-

(^ XiE^^^) + ^{zia - 1 )-^ ^

2=1

2=1

[e^

1 ]

From a nonlinear system of equations U{y;@) = 0, we obtain the iterative procedure of the
EM-algorithm as

a

(i+i) -

—n

nC'(0(j+i)) ^

n

—iY

/5(i+i) yO')

[^a+i)] 2 '^

'v(j) 7’ •

- n

) + E(

z,-a

2=1

[e

= 0 ,

- 1 ]

nx +

^^ 0 + 1 ) 3 ;.

— n) —

Y'>

70 + 1 )

' i=l

(Y^ Xie^'

0 + 1 ).

+E(

ZM''-

(j) _

_ I'j I ^

+ 1 + 1)12 J-J W

vO + l)

2=1

lp+ 0 ' + l)a;i 11

^- 1 ]

= 0 ,

where ^(i+i) and yb’+i) are found numerically. Here, for i = 1, 2,... , n, we have that

Za — 1 “b

where [1 — e

We can use the results of

11 ,-s

Louii

(jl982l l to obtain the standard errors of the estimators from
the EM-algorithm. Consider lc{G;x) = E{Ic{&;y)\x), where Ic{Q;y) = is the

4x4 observed information matrix.If lm{®] x) = Var[U{y, 0)|a;], then, we obtain the observed
information as

/(©; x) = lc{&-, x) - lm{®; x).

The standard errors of the MLEs of the EM-algorithm are the square root of the diagonal
elements of the/(0; a;). The computation of the se matrices are too long and t edious. Therefore,

we did not present the details. Reader can see
these values.

Mahmoudi and Jafaril (120121 ) how to calculate

6 Simulation study

We performed a simulation in order to investigate the proposed estimator of a, /3, 7 and 9 of the
proposed EM-scheme. We generated 1000 samples of size n from the GGG distribution with
/3 = 1 and 7 = 0.1. Then, the averages of estimators (AE), standard error of estimators (SEE),
and averages of standard errors (ASE) of MLEs of the EM-algorithm determined though the
Fisher information matrix are calculated. The results are given in Table [21 We can find that

(i) convergence has been achieved in all cases and this emphasizes the numerical stability of
the EM-algorithm,

(ii) the differences between the average estimates and the true values are almost small,

(hi) the standard errors of the MLEs decrease when the sample size increases.

Table 2; The average MLEs, standard error of estimators and averages of standard errors for

the GGG distribution.

parameter

AE

SEE

ASE

n

a

e

a

/?

7

e

a

/?

7

e

a

/?

7

e

50

0.5

0.2

0.491

0.961

0.149

0.204

0.114

0.338

0.265

0.195

0.173

0.731

0.437

0.782

0.5

0.5

0.540

0.831

0.182

0.389

0.160

0.337

0.260

0.263

0.210

0.689

0.421

0.817

0.5

0.8

0.652

0.735

0.154

0.684

0.304

0.377

0.273

0.335

0.309

0.671

0.422

0.896

1.0

0.2

0.988

0.972

0.129

0.206

0.275

0.319

0.191

0.209

0.356

0.925

0.436

0.939

1.0

0.5

1.027

0.852

0.147

0.402

0.345

0.352

0.226

0.283

0.408

0.873

0.430

0.902

1.0

0.8

1.210

0.711

0.178

0.745

0.553

0.365

0.230

0.342

0.568

0.799

0.433

0.898

2.0

0.2

1.969

0.990

0.084

0.216

0.545

0.305

0.151

0.228

0.766

1.135

0.422

0.902

2.0

0.5

1.957

0.842

0.113

0.487

0.608

0.334

0.192

0.277

0.820

1.061

0.431

0.963

2.0

0.8

2.024

0.713

0.161

0.756

0.715

0.396

0.202

0.353

1.143

0.873

0.402

0.973

100

0.5

0.2

0.491

0.977

0.081

0.212

0.084

0.252

0.171

0.179

0.125

0.514

0.283

0.561

0.5

0.5

0.528

0.883

0.109

0.549

0.124

0.275

0.178

0.247

0.155

0.504

0.275

0.567

0.5

0.8

0.602

0.793

0.136

0.769

0.215

0.323

0.194

0.299

0.220

0.466

0.259

0.522

1.0

0.2

0.974

0.997

0.102

0.226

0.195

0.242

0.129

0.206

0.251

0.645

0.280

0.767

1.0

0.5

1.030

0.875

0.113

0.517

0.262

0.291

0.155

0.270

0.298

0.651

0.295

0.843

1.0

0.8

1.113

0.899

0.117

0.846

0.412

0.342

0.177

0.331

0.400

0.600

0.287

0.781

2.0

0.2

1.952

0.995

0.138

0.221

0.424

0.237

0.117

0.209

0.524

0.922

0.321

0.992

2.0

0.5

2.004

0.885

0.110

0.518

0.493

0.283

0.131

0.274

0.601

0.873

0.321

0.966

2.0

0.8

2.028

0.981

0.104

0.819

0.605

0.350

0.155

0.339

0.816

0.717

0.289

0.946

7 Real examples

In this Section, we consider two real data sets and fit the Gompertz, GGG, GGP, GGB (with
m = 5), and GGL distributions. The first data set is negatively skewed, and the second data set
is positively skewed, and we show that the proposed distributions fit both positively skewed and
negatively skewed data well. Eor each data, the MLE of parameters (with standard deviations)
for the distributions are obtained. To test the goodness-of-fit of the distributions, we calculated
the maximized log-likelihood, the Kolmogorov-Smirnov (K-S) statistic with its respective p-
value, the AIG (Akaike Information Griterion), AIGG (AIG with correction), BIG (Bayesian
Information Criterion), CM (Cramer-von Mises statistic) and AD (Anderson-Darling statistic)
for the six distributions. Here, the significance level is 0.10. To show that the likelihood

Table 3: The strengths of glass fibers.

0.55,

0.93,

1.25,

1.36,

1.49, 1.52,

1.58,

1.61, 1.64,

1.68, 1.73,

1.81, 2 . 00 ,

0.74,

1.04, 1.27,

1.39,

1.49,

1.53,

1.59,

1.61, 1 . 66 ,

1 . 68 ,

1.76, 1.82,

2.01, 0.77,

1 . 11 , 1.28,

1.42,

1.50, 1.54,

1.60,

1.62,

1 . 66 ,

1.69,

1.76, 1.84,

2.24,

0.81, 1.13,

1.29, 1.48,

1.50, 1.55,

1.61,

1.62, 1 . 66 ,

1.70,

1.77,

1.84,

0.84,

1.24, 1.30,

1.48,

1.51, 1.55,

1.61, 1.63,

1.67, 1.70,

1.78,

1.89

Table 4: The phosphorus concentration in

the leaves.

0 . 22 ,

0.17,

0 . 11 ,

0 . 10 ,

0.15, 0.06,

0.05,

0.07, 0.12,

0.09, 0.23,

0.25, 0.23,

0.24,

0 . 20 , 0.08

0 . 11 ,

0 . 12 ,

0 . 10 ,

0.06,

0.20, 0.17,

0 . 20 ,

0 . 11 , 0.16,

0.09, 0.10,

0 . 12 , 0 . 12 ,

0 . 10 ,

0.09, 0.17

0.19,

0 . 21 ,

0.18,

0.26,

0.19, 0.17,

0.18,

0.20, 0.24,

0.19, 0.21,

0.22, 0.17,

0.08,

0.08, 0.06

0.09,

0 . 22 ,

0.23,

0 . 22 ,

0.19, 0.27,

0.16,

0.28, 0 . 11 ,

0 . 10 , 0 . 20 ,

0.12, 0.15,

0.08,

0.12, 0.09

0.14,

0.07,

0.09,

0.05,

0.06, 0 . 11 ,

0.16,

0.20, 0.25,

0.16, 0.13,

0 . 11 , 0 . 11 ,

0 . 11 ,

0.08, 0.22

0 . 11 ,

0.13,

0 . 12 ,

0.15,

0 . 12 , 0 . 11 ,

0 . 11 ,

0.15, 0.10,

0.15, 0.17,

0.14, 0.12,

0.18,

0.14, 0.18

0.13,

0 . 12 ,

0.14,

0.09,

0.10, 0.13,

0.09,

0 . 11 , 0 . 11 ,

0.14, 0.07,

0.07, 0.19,

0.17,

0.18, 0.16

0.19,

0.15,

0.07,

0.09,

0.17, 0.10,

0.08,

0.15, 0.21,

0.16, 0.08,

0 . 10 , 0.06,

0.08,

0.12, 0.13

equations have a unique solution in the parameters, we plot the profile log-likelihood functions
of /3, 7 , a and Q for the six distributions.

First, we con s ider t he data consisting of the strengths of 1.5 cm glass fibers given in

Smith and Navlon (llQSTIi and measured at t

data is also studied by

le b

Barreto-Souza et al.

ational Physical Laboratory, England. This

( 201C ) and is given in Table [3l

The results are given in Table [5] and show that the GGG distribution yields the best fit
among the GGP, GGB, GGL, GG and Gompertz distributions. Also, the GGG, GGP, and
GGB distribution are better than GG distribution. The plots of the pdfs (together with the
data histogram) and cdfs in Figure [7] confirm this conclusion. Figures [9] show the profile log-

likelihood functions of /3, 7 , a and d for the six distrib utions. _

As a second example, we consider a data set from iFonseca and Francal (j2007l i. who stud¬

ied the soil fertility influence and the characterization of the biologic fixation of N 2 for the
Dimorphandra wilsonii rizz growth. For 128 plants, they made meas u res o f the phosphorus
concentration in the leaves. This data is also studied by

Silva et al.

(l2013l i and is given in

Table 01 Figures [TO] show the profile log-likelihood functions of (3, 7 , a and 9 for the six
distributions.

The results are given in Table [ 6 l Since the estimation of parameter 9 for GGP, GGB,
and GGL is close to zero, the estimations of parameters for these distributions are equal to
the estimations of parameters for GG distribution. In fact. The limiting distribution of GGPS
when 0 —)• 0^ is a GG distribution (see Proposition [2|) . Therefore, the value of maximized log-
likelihood, log(L), are equal for these four distributions. The plots of the pdfs (together with

the data histogram) and cdfs in Figure [8] confirm these conclusions. Note that the estimations
of parameters for GGG distribution are not equal to the estimations of parameters for GG
distribution. But the log(L)’s are equal for these distributions. However, from Table [6] also we
can conclude that the GG distribution is simpler than other distribution because it has three
parameter but GGG, GGP, GGB, and GGL have four parameter. Note that GG is a special
case of GGPS family.

Histogram

Empirical Distribution

X

Figure 7: Plots (pdf and cdf) of fitted Gompertz, generalized Gompertz, GGG, GGP, GGB
and GGL distributions for the first data set.

Histogram

Empirical Distribution

X

Figure 8: Plots (pdf and cdf) of fitted Gompertz, generalized Gompertz, GGG, GGP, GGB
and GGL distributions for the second data set.

Table 5; Parameter estimates (with std.), K-S statistic, p-value, AIC, AICC and BIC for the
hrst data set,_

Distribution

Gompertz

GG

GGG

GGP

GGB

GGL

/3

0.0088

0.0356

0.7320

0.1404

0.1032

0.1705

s.e.(/3)

0.0043

0.0402

0.2484

0.1368

0.1039

0.2571

7

3.6474

2.8834

1.3499

2.1928

2.3489

2.1502

s.e.{j)

0.2992

0.6346

0.3290

0.5867

0.6010

0.7667

a

—

1.6059

2.1853

1.6205

1.5999

2.2177

s.e.{a)

—

0.6540

1.2470

0.9998

0.9081

1.3905

e

—

—

0.9546

2.6078

0.6558

0.8890

s.e.{6)

—

—

0.0556

1.6313

0.5689

0.2467

-log(L)

14.8081

14.1452

12.0529

13.0486

13.2670

13.6398

K-S

0.1268

0.1318

0.0993

0.1131

0.1167

0.1353

p-value

0.2636

0.2239

0.5629

0.3961

0.3570

0.1992

AIC

33.6162

34.2904

32.1059

34.0971

34.5340

35.2796

AICC

33.8162

34.6972

32.7956

34.78678

35.2236

35.9692

BIC

37.9025

40.7198

40.6784

42.6696

43.1065

43.8521

CM

0.1616

0.1564

0.0792

0.1088

0.1172

0.1542

AD

0.9062

0.8864

0.5103

0.6605

0.7012

0.8331

Table 6; Parameter estimates (with std.), K-S statistic, p-value, AIC, AICC and BIC for the
second data set._

Distribution

Gompertz

GG

GGG

GGP

GGB

GGL

/3

1.3231

13.3618

10.8956

13.3618

13.3618

13.3618

s.e.0)

0.2797

4.5733

8.4255

5.8585

6.3389

7.3125

7

3.1500

4.0158

3.1500

3.1500

3.1500

s.e.{j)

2.1865

3.6448

2.4884

2.6095

2.5024

a

—

6.0906

5.4236

6.0906

6.0906

6.0905

s.e.(d)

—

2.4312

2.8804

2.6246

2.7055

2.8251

e

—

—

-0.3429

1.0 X 10"®

1.0 X 10"®

1.0 X 10"®

s.e.{6)

—

—

1.2797

0.8151

0.2441

0.6333

- log(T)

-184.5971

-197.1326

-197.1811

-197.1326

-197.1326

-197.1326

K-S

0.1169

0.0923

0.0898

0.0923

0.0923

0.0923

p-value

0.06022

0.2259

0.2523

0.2259

0.2259

0.2259

AIG

-365.1943

-388.2653

-386.3623

-386.2653

-386.2653

-386.2653

AIGG

-365.0983

-388.0717

-386.0371

-385.9401

-385.9401

-385.9401

BIG

-359.4902

-379.7092

-374.9542

-374.8571

-374.8571

-374.8571

CM

0.3343

0.1379

0.1356

0.1379

0.1379

0.1379

AD

2.3291

0.7730

0.7646

0.7730

0.7730

0.7730

Log Likelihood Lqq Likelihood

Figure 9: The profile log-likelihood functions for Gompertz, generalized Gompertz, GGG, GGP,
GGB and GGL distributions for the first data set.

Log Likelihood Lqq Likelihood

Figure 10: The profile log-likelihood functions for Gompertz, generalized Gompertz, GGG,

GGP, GGB and GGL distributions for the second data set.

Appendix

We demonstrate those parameter intervals for which the hazard function is decreasing, i ncreas¬
i ng an d bathtub shaped, and in order to do so, we follow closely a theorem given by

Glaser

(j 19801 1. Define the function t{x) = where f'{x) denotes the first derivative of f{x) in

(12.3p . To simplify, we consider u = 1 — exp(^(e'>'^ “ !))•

A.l

Consider the GGG hazard function in (14.2h . then we define

/ \ 1 — a 2a9u°^~^

"W = 7(T = — +

If a > 1, then t '{ u ) > 0, and h{.) is an increasing function. If 0 < a < 1, then

, , 2a6'^ 1

lim r (u) = -oo, lim r (u) = + (a - 1)(1 - > 0.

Since the limits have different signs, the equation t '{ u ) = 0 has at least one root. Also, we
can show that t "{ u ) > 0. Therefore, the equation t '{ u ) = 0 has one root. Thus the hazard
function is decreasing and bathtub shaped in this case.

A.2

The GGP hazard rate is given by h{u) = 6a(5u'^ — e®““). We define = log[/i(n)].

Then, its hrst derivative is

/ / X a — I ^ n M

Tj {u) =-h aOe^

OL—1

U

g6» _

It is clearly for a > 1, r]'{u) > 0 and h{u) is increasing function. If 0 < a < 1, then

lim rj'{u) = —oo, lim rj'{u) = 0,

So the equation t'{u) = 0 has at least one root. Also, we can show that r"(u) > 0. It implies
that equation r]'(u) = 0 has a one root and the hazard rate increase and bathtub shaped.

B.

B.l

Let to,(a) = Eti = W Eti log(C'(«i“)). For GGG,

wi{a) = 26

2=1

tf log tj dwi{a)

1 — 9tf ’ da

» Li _ g^al

2=1

For GGP,

For GGL,

rt;i(a) = O'^tflogU,

i=l

wi{a) = 9

dwi{a)
da

tf log ti dwi{a)

> 0 .

2=1

2=1

1 - 6»i“ ’

da

2=1

1 - 9t^

For GGB,

u,(a) = (m - 1)9 V = {.» - 1)9 E ‘f [^1" > 0.

i=l

2=1

Therefore, rt;i(a) is strictly increasing in a and

lim gi{a;/3,-f,9,x) = oo, lim gi(a;/?, 7 ,0, x) = y^logiU).

Q->- 0 + a^oo

1=1

Also,

n ^ ^ ji 6C"(6) ^ ^

gi(a;^,7,6',x) < - + ^log{ti), gi{a;/ 3 ,j,e,x) > - + ( + l)^log{ti).

“ i=i a o

Hence, gi(a; 7 , 0, x) < 0 when ^ + 1]?=! log(^*) < 0i and gi(a;/3, 7 , 6 », x) > 0 when | +
Sr=i log(^*) > 0 - The proof is completed.

B.2

It can be easily shown that

-1 ”

lim g 2 (/ 3 ;a, 7 , 6 »,x) = 00 , lim g 2 (/ 3 ; a, 7 , 0 , x) = —- 1).

/3-)^0+ /9-)^oo 7 ^

^=l

Since the limits have different signs, the equation g 2 (/ 3 ; a, 7 ,0, x) = 0 has at least one root with
respect to /3 for fixed values a, 7 and 9. The proof is completed.

B.3

a) For GGP, it is clear that

limg 3 ( 6 »;a,/ 3 , 7 ,x) = lim g 3 ( 6 »; a,/?, 7 , x) = - 00 .

e-l-O ^^ I d^oo

i=l

Therefore, the equation g 3 ( 0 ; a, /3, 7 , x) = 0 has at least one root for 0 > 0, if Yll=i “ f > 0

or > §•

b) For GGG, it is clear that

lim g 3 ( 6 »;a,/ 3 , 7 ,x) = - 00 ,
0^00

n

lim g 3 (( 9 ; 0 ,^, 7 , x)
0-)-O+

= -n + 2Y,tt

2 = 1

Therefore, the equation g3{9, (3, 7 , x) = 0 has at least one root for 0 < 0 < 1, if —n+2 >

0 or > §•

For GGL, it is clear that

liing 3 ( 6 »;a,/ 3 , 7 ,x) = Vtf - lim g 3 ( 6 '; a,/3, 7 , x) = - 00 .

2=1

2’ 0^1 ‘

Therefore, the equation g 3 ( 0 ; a, (3, 7 , x) = 0 has at least one root for 0 < 0 < 1, if ^11=1 “ § >

Oor >§•

For GGB, it is clear that

liing 3 (p;a,/ 3 , 7 ,x) = 1) - , lini g 3 (p; a, 7 , x) ^ ^ 1

p—>-0 ^^ z p —>0 ^

i=l i=l

Therefore, the equation g 3 (p; a, /3, 7 , x) = 0 has at least one root for 0 < p < 1, if XlILi tfim —

m-

~T

_ n(m-l)

> 0 and E?.i < 0 or E?.i <f > ? and E?=i *r" > iSk-

c.

Consider t* = 1 — e 7 ^*^^ ' Then, the elements of 4 x 4 observed information matrix In{Q)
are given by

laa or 9 —

oa‘^ a

n/^ < ^ I

2=1

Iaj3 —

d'^L

dad 13 ^ ■ 7

"■ g7d;i _ 2 0 ”■

E

+ -111(0 iog(t.) +1)§^

^ 2=1 ^ ^

^ 0:7 —

+Q;6»tf log(ti)

d^ln
dadj

C'"{9tf)C'{etf) - {C''{9tf)f

{C'{9tf)f

e7dir(7x. _i) + i ep

2=1

7"

] + 4 + l][(alog(t,) + 1 )

3-a9

I(3f} =

Ip'y —

+a9ti log(ti)

dX

dad9

C"'{9tf)C'{9tf) - {C''{9tf))

^ ^=l

ll(Q+a\\2

C''{9tt)

C'{9tf)

2=1

ic'rnw

C"{9tf) C'''{9tf)C’{9tf) - {C"{9t<t)f

C'{9tf) *

{C'{9tf)r

—n

dX _

dl3^ /32

dX {a - 2 )

+ 9a‘^'^tj

^^^a^tf)C'{9tf) - iC"{9tf)y

2=1

n

7 ^ X'{9tf)

{C'{9tf)Y

dpd'y

E(*""-(7a^. -D + D + aOY. - ') + ‘>1?^

2=1

^ C'''i9tf)C'i9X-iC''i9tf)y

X ^ X{9tf)W

Ifie =

-^77 ~

lO'y

lee

dHn

d/3de

E*:

2 = 1

l(aj.a\'\2

ic'ietf))

rfil OR JE —2 q.'^p'Y^i

- 1) + 1] + (a - - 1) + 1) + ^1

2=1

-2

2=1

■

/-yO

7 metf)

7 '

+ ^(e^"*(7x.-l) + l)
T

dHn
d9dj

C"(0tf)

( C "( 02 f )))2

_ 1) +1,1^ + - (a”(«.“))^

i=l ^

C"(0tf

502 02

-n ^.2a. C'"{9tf)C'{0tf)-{C"{0t<l)r
ft2 Z-,- (riltQ+a\\2 J

2=1

(C'(0tf))2

(C'(02f))2

C''{0)C'{0) - (C'(0))2
(C"(0))2 ^

Acknowledgements

The authors would like to thank the anonymous referees for many helpful comments and
suggestions.

References

Adamidis, K., Dimitrakopoulou, T., and Loukas, S. (2005). On an extension of the exponential-
geometric distribution. Statistics and Probability Letters, 73(3);259-269.

Adamidis, K. and Loukas, S. (1998). A lifetime distribution with decreasing failure rate.
Statistics and Probability Letters, 39(l):35-42.

Barreto-Souza, W., Morals, A. L., and Cordeiro, G. M. (2011). The Weibull-geometric distri¬
bution. Journal of Statistical Computation and Simulation, 81(5):645-657.

Barreto-Souza, W., Santos, A. H. S., and Cordeiro, G. M. (2010). The beta generalized expo¬
nential distribution. Journal of Statistical Computation and Simulation, 80(2):159-172.

Cancho, V. G., Louzada-Neto, F., and Barriga, G. D. C. (2011). The Poisson-exponential
lifetime distribution. Computational Statistics and Data Analysis, 55(l):677-686.

Casella, G. and Berger, R. (2001). Statistical Inference. Duxbury, Pacific Grove, California,

USA.

Chahkandi, M. and Ganjali, M. (2009). On some lifetime distributions with decreasing failure
rate. Computational Statistics and Data Analysis, 53(12) :4433-4440.

Dempster, A. P., Laird, N. M., and Rubin, D. B. (1977). Maximum likelihood from incomplete
data via the EM algorithm. Journal of the Royal Statistical Society. Series B (Methodologi¬
cal), 39(l):l-38.

El-Gohary, A., Alshamrani, A., and Al-Otaibi, A. N. (2013). The generalized Gompertz distri¬
bution. Applied Mathematical Modelling, 37(l-2):13-24.

Elores, J., Borges, P., Cancho, V. G., and Louzada, F. (2013). The complementary exponential
power series distribution. Brazilian Journal of Probability and Statistics, 27(4);565-584.

Fonseca, M. and Franca, M. (2007). A influencia da fertilidade do solo e caracterizagao da
fixagao biologica de N 2 para o crescimento de. Dimorphandra wilsonii rizz. Masters thesis,
Universidade Federal de Minas Gerais.

Glaser, R. E. (1980). Bathtub and related failure rate characterizations. Journal of the Amer¬
ican Statistical Association, 75(371);667-672.

Gupta, R. D. and Kundu, D. (1999). Generalized exponential distributions. Australian & New
Zealand Journal of Statistics, 41(2):173-188.

Ku§, C. (2007). A new lifetime distribution. Computational Statistics and Data Analysis,
51(9):4497-4509.

Lawless, J. F. (2003). Statistical Models and Methods for Lifetime Data. Wiley-Interscience,
second edition.

Louis, T. A. (1982). Finding the observed information matrix when using the EM algorithm.
Journal of the Royal Statistical Society. Series B (Methodological), 44(2):226-233.

Louzada, F., Roman, M., and Cancho, V. G. (2011). The complementary exponential geometric
distribution; Model, properties, and a comparison with its counterpart. Computational
Statistics and Data Analysis, 55(8);2516-2524.

Louzada-Neto, F., Cancho, V. G., and Barriga, G. D. C. (2011). The Poisson-exponential
distribution; a Bayesian approach. Journal of Applied Statistics, 38(6);1239-1248.

Mahmoudi, E. and Jafari, A. A. (2012). Generalized exponential-power series distributions.
Computational Statistics and Data Analysis, 56(12) :4047-4066.

Mahmoudi, E. and Jafari, A. A. (2014). The compound class of linear failure rate-power series
distributions: model, properties and applications. arXiv preprint arXiv:1402.5282.

Marshall, A. W. and Olkin, I. (1997). A new method for adding a parameter to a fam¬
ily of distributions with application to the exponential and Weibull families. Biometrika,
84(3):641-652.

Morais, A. L. and Barreto-Souza, W. (2011). A compound class of Weibull and power series
distributions. Computational Statistics and Data Analysis, 55(3):1410-1425.

Noack, A. (1950). A class of random variables with discrete distributions. The Annals of
Mathematical Statistics, 21(1):127-132.

Shannon, C. (1948). A mathematical theory of communication. Bell System Teehnical Journal,
27:379-432.

Silva, R. B., Bourguignon, M., Dias, C. R. B., and Cordeiro, G. M. (2013). The compound class
of extended Weibull power series distributions. Computational Statisties and Data Analysis,
58:352-367.

Smith, R. L. and Naylor, J. G. (1987). A comparison of maximum likelihood and Bayesian
estimators for the three-parameter Weibull distribution. Applied Statisties, 36(3):358-369.

Tahmasbi, R. and Rezaei, S. (2008). A two-parameter lifetime distribution with decreasing
failure rate. Computational Statistics and Data Analysis, 52(8):3889-3901.

Hazard

0 12 3 4

X

Hazard

0 12 3 4

X

Hazard

0 12 3 4

X

Hazard

0 12 3 4

X